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- DISCRETE TIME SIGNALS 


A number of important types of signals and their operations are introduced. 

Linear and shift-invariant systems are discussed mostly because they are easier to analyze and 
implement. The convolution and the difference equation representations are given special attention 
because of their importance in digital signal processing and in MATLAB. 


Signals are broadly classified into analog and discrete signals. An analog signal will be denoted by 
x a (t), in which the variable t can represent any physical quantity, but we will assume that it 
represents time in seconds. A discrete signal will be denoted by x(n), in which the variable n is 
integer- valued and represents discrete instances in time. Therefore it is also called a discrete-time 
signal, which is a number sequence and will be denoted by one of the following notations: 

x(n) = {x(n)} = {. . . , x( - 1 ), x(0), x{1), . . .} 

t 

Where the up-arrow indicates the sample at n =0. 

A correct representation of x(n) would require two vectors, one each for x and n. For example, a 
sequence x(n) = {2., 1,-1, 0, 1, 4, 3, 7} can be represented in MATLAB by: 

t 

» n=[-3,-2,-1,0,1,2,3,4]; x=[2,1 ,-1,0,1,4,3,7]; 


TYPES OF SEQUENCES 

There are several elementary sequences in digital signal processing for analysis purposes. Their 
definitions and MATLAB representations follow: 


1. Unit sample sequence 




= <_ 0 , 0 , 1 , 0 , 0 ,... 


In MATLAB the function zeros(1 ,N) generates a row vector of N zeros, which can be used to 
implement 6(n) over a finite interval. However, the logical relation n==0 is an elegant way of 
implementing 6 (n). For example, to implement 



n = no 
n 4 no 


Over the m <m <m interval, we will use the following MATLAB function: 


function [x,n] = impseq(nO, nl, n2) 

% Generates x(n) = delta(n-nO); nl 

o 

o 


% [x,n] = impseq(nO , nl,n2) 

o 

o 


== 0 ]; 


<= n <= n2 


n 


[nl:n2]; x 


[(n-nO) 



2. Unit step sequence 



1 > n > 0 
0 . n < 


1 


In MATLAB the function ones(1 ,N) generates a row vector of N ones. It can be used to generate u(n) 
over a finite interval. Once again an elegant approach is to use the logical relation n>=0. To 
implement: 


, , ft n>n o 

li(n-no) =-f n 

I U s n < no 


Over the m <m <m interval, we will use the following MATLAB function. 

function [x,n] = stepseq(nO,nl,n2) 

% Generates x(n) = u(n-nO); nl <= n <= n2 

o 

o 

% [x,n] = stepseq(nO , nl,n2) 

o 

o 

n = [nl:n2]; x = [(n-nO) >= 0]; 


3. Real-valued exponential sequence 

x(n) = a", Vn\ a e R 

In MATLAB an array operator “, A ” is required to implement a real exponential sequence. For example, 
to generate x(n) = (0. 9 )n, 0 <n <10, we will need the following MATLAB script: 

>>n= [0:10]; x= (0.9). A n; 

4. Complex-valued exponential sequence 

x (n) = d a+jm) n Vn 

Where a produces an attenuation (if <0) or amplification (if 20) and coo is the frequency in radians. A 
MATLAB function exp is used to generate exponential sequences. For example, to generate 

4tf)=exp[(2+ j3)n], 0, we will need the following MATLAB script: 

>> n = [0:10]; x = exp ((2+3j)*n); 

5 . Sinusoidal sequence 


x(n) = A COS ( co 0 n + 0 a ), Vn 



Where A is an amplitude and is the phase in radians. A MATLAB function cos (or sin) is used to 
generate sinusoidal sequences. 


For example, to generate x(ri) = 3cos(0.1 an + n/3) + 2 sin(0. 5 nri), 
0 <n <10, we will need the following MATLAB script: 


» n = [0:10]; x = 3*cos(0.1*pi*n+pi/3) + 2*sin(0.5*pi*n); 

6. Random sequences 

Many practical sequences cannot be described by mathematical expressions like those above. These 
sequences are called random (or stochastic) sequences and are characterized by parameters of the 
associated probability density functions. In MATLAB two types of (pseudo-) random sequences are 
available. The rand(1,N) generates a length N random sequence whose elements are uniformly 
distributed between [0, 1]. The randn(1,N) generates a length N Gaussian random sequence with 
mean 0 and variance 1. Other random sequences can be generated using transformations of the 
above functions. 

7. Periodic sequences 

A sequence x(n) is periodic if x(n) = x(n + A), Vn. The smallest integer N that satisfies this relation is 
called the fundamental period. We will use ~x(ri) to denote a periodic sequence. 

To generate P periods of ~x(n) from one period {x(n), 0 <n <N-\ }, we can copy x{n) P times: 


» xtilde = [x,x,...,x]; 

But an elegant approach is to use MATLAB’s powerful indexing capabilities. 

First we generate a matrix containing P rows of x(n) values. 

Then we can concatenate P rows into a long row vector using the construct (:). However, this 
construct works only on columns. Hence we will have to use the matrix transposition operator ’ to 
provide the same effect on rows. 

>> xtilde = x' * ones(l,P); % P columns of x; x is a row vector 
>> xtilde = xtilde(:); % long column vector 

>> xtilde = xtilde'; % long row vector 


- Operations on Signals 

1. Signal addition 

The sample-by-sample addition is given by: 

{Mn)}+ {x 2 (ri)}= fa (n) + x 2 (n)} 

The lengths of xi(n) and xz{n) must be the same. If sequences are of unequal lengths, or if the sample 
positions are different for equal length sequences, then we cannot directly use the operator +. 


We need to augment X\ (n) and x 2 (ri) in order to make them at the same position vector n ( hence 
the same length. This requires careful attention to the Matlab indexing operations. In particular, 
logical operation of intersection relational operations like and and the find function are 
required to make x\ (n) and xi(n) of equal length. The following function, called the sigadd function, 
demonstrates these operations. 

function [y,n] = sigadd(xl,nl,x2,n2) 

% implements y(n) = xl(n)+x2 (n) 

o 

o 

% [y,n] = sigadd(xl, nl,x2, n2) 

% y = sum seguence over n, which includes nl and n2 
% xl = first seguence over nl 

% x2 = second seguence over n2 (n2 can be different from nl) 

o 

o 

n = min(min(nl),min(n2)):max(max(nl),max(n2)); % duration of y(n) 
yl = zeros(1,length(n)); y2 = yl; % initialization 

yl(find((n>=min(nl))&(n<=max(nl))==1))=xl; % xl with duration of y 
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2; % x2 with duration of y 
y = yl+y2; % seguence addition 


Example 1: 

Let Xi(n) = {1,2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2,1}. 

r 

And x 2 (n) = {3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3}. 

r 


Determine and plot the following sequences 
x 3 (n) = x ; (n) + x 2 {n) 


Solution 1: 

clear, clc, elf 

xl = [1, 2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2,1] ; nl = -5 : 7; 

x2= [ 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3 ]; n2 = -l:9; 

[y,n] = sigadd(xl,nl,x2,n2); 
subplot (2,2,1) 

stem(nl,xl)/title( ’signal xl 1 );ylabel( 1 xl(nl) ! );xlabel( ’nl’ )/axis([-5 7 0 

8 ] ) ; 

subplot (2,2,2) 

stem(n2,x2) /title ( ’signal x2’ )/ylabel( 1 x2(n2 ) ’ )/xlabel( f n2')/axis([-l 9 0 

8 ] ) / 

subplot (2,2,3:4) 

stem(n,y)/title ( ’signal y' )/ylabel( 'y(n)' )/xlabel( ! n’ )/axis([-5 9 0 11])/ 
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figure 1.1 : Solution Results of Example 1 


2. Signal multiplication 

This is a sample-by-sample (or “dot”) multiplication) given by 


{xi(n)} {x 2 (n)} = {xi(n) x 2 (n)} 

It is implemented in MATLAB by the array operator Once again, the similar restrictions apply for the 
.* operator as for the + operator. Therefore we have developed the sigmult function, which is similar to 
the sigadd function. 

Exercise 1: 

Write a function in Matlab named (( sigmul.m )) that performs the function of signal multiplication. 


Solution el: 

function [y,n] = sigmult(xl,nl,x2,n2) 

% implements y(n) = xl(n)*x2(n) 

o 

o 

% [y,n] = sigmult(xl,nl,x2 , n2 ) 

% y = product sequence over n, which includes nl and n2 
% xl = first sequence over nl 

% x2 = second sequence over n2 (n2 can be different from nl) 

o 

o 

n = min(min(nl) , min(n2)):max(max(nl) , max(n2)); % duration of y(n) 
yl = zeros(1,length(n)); y2 = yl; % 

yl(find((n>=min(nl))&(n<=max(nl))==1))=xl; % xl with duration of y 
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2; % x2 with duration of y 
y = yl .* y2; % sequence multiplication 

















































3. Signal scaling 

In this operation each sample is multiplied by a scalar a. 

a fin)} = {a 4n)} 

An arithmetic operator (*) is used to implement the scaling operation in MATLAB. 

4. Signal shifting 

In this operation, each sample of x(n) is shifted by an amount k to obtain a shifted sequence y(n). 

y(n) = {x(n - k)} 

If we let m = n-k, then n = m+k and the above operation is given by 

y{m + k) = {x (m)} 

Hence this operation has no effect on the vector x, but the vector n is changed by adding k to each 
element. This is shown in the function sigshift. 

function [y,n] = sigshift(x,m,k) 

% implements y(n) = x(n-k) 

o 

o 

% [y,n] = sigshift(x,m,k) 

o 

o 

n = m+k; y = x; 


Example 2: 


Let x(n) = {1,2, 3, 4, 5}. 

t 

Find and plot the following: 


a- xl (n)= x(n-2) 
b- x2 (n)= x(n+3) 

Solution 2: 


% delay with 2 sec 
% advance with 3 sec 


clear all , clc 
x = [ 1 2 3 4 5 ] ; 
m = [—1:3]; 
kl = 2; 
k2 = -3; 


% Original Signal 
% Original position 
% value of the delay for signal 1 
% Value for the delay for signal 2 


[xl,nl] = sigshift(x,m,kl); 
[x2,n2] = sigshift(x A m A k2); 


Subplot ( 3,1,1) 

stem(m,x); title ( Original signal x'); xlabel(’m ? ); ylabel(' x(m) ' ); axis 
([-5 5 0 8]) ; 


Subplot ( 3,1,2) 

Stem(nl,xl); title ( f Delayed signal xl ! ) ; xlabel ( ’ nl ' ) ; ylabel ( xl (nl) ’ ) ; 
axis ([-5 5 0 8]) ; 


Subplot ( 3,1,3) 

Stem(n2,x2); title ( Advanced signal x2' ) ; xlabel ( ' n2 ! ) ; 
ylabel( 'x2(n2)' );axis ([-5 5 0 8]); 
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figure 1.1 : Solution results of Example 2 


5. Signal folding 

In this operation each sample of x(n) is flipped around n = 0 to obtain a folded sequence y{n). 

y{n) = {x( -n)} 

In MATLAB this operation is implemented by fliplr(x)function for sample values and by -fliplr(n) 
function for sample positions as shown in the sigfold function: 

function [y,n] = sigfold(x,n) 

% implements y(n) = x(-n) 

o 

o 

% [y,n] = sigfold(x,n) 

o 

o 

y = fliplr(x); n = -fliplr(n); 


6. Sample summation 


This operation differs from signal addition operation. It adds all sample values of x(n) between m and 

HI. 

n 2 

v 

7 x(n ) = x(ni) + * * * + x(n 2 ) 

n- n x 


It is implemented by the sum(x (nl :n2)) function. 

















7. Sample product 


This operation also differs from signal multiplication operation. It multiplies all sample values of x{n) 
between m and m. 


n 2 

T I x(n) = x(m) + • • • + x(n2) 

«1 


It is implemented by the prod(x(n1:n2)) function. 

8. Signal energy 

The energy of a sequence x{n) is given by 



DO DO 

:E(ti);r*(n.) = |:r{n) 

—oo —oo 



2 


Where superscript ^denotes the operation of complex conjugation. 
The energy of a finite-duration sequence x(ri) can be computed in 
MATLAB using 

>> Ex = sum(x .* conj (x)); % one approach 
>> Ex = sum(abs(x) . A 2); % another approach 


9. Sample power 

The average power of a periodic sequence ~x(n) with fundamental period N is given by: 




9 





Exercise 2: 


Generate and plot each of the following sequences over the indicated interval. 

a) x(n) = 25(n + 2) - 5(n - 4), -5 < n < b. 

b) x(n) = n[u(n)-u(n-^0)]+^Oe-o.3(n-^o)[u(n-^0)-u(n-20)], 0 < n < 20. 

c) x(n) = cos(0.04nn) + 0.2 w(n), 0 < n < 50, 



where w(n) is a Gaussian 

random sequence with zero mean and unit variance. 
~x(n) = 5, 4, 3, 2, 1, 5, 4, 3, 2, 1,5, 4, 3, 2, 1, ... 7 ; 

t 


-10 < n <9. 


Solution e2: 

a. x(ri) = 25(n + 2) - 5(n -4), -5 < n < 5. 

n = [—5:5]; 

x = 2*impseq (-2, -5,5) - impseq(4,-5,5); 

subplot (2,2,1); stem(n,x)/ title(’ Sequence in Exercise a') 
xlabel( 'n' ); ylabel( ! x(n)’ ); 

b. x(n) = n[u(n)-u(n^O)]+'\Oe-o.3(n- j \o)[u(ri^O)-u(n-20)], 0 < n < 20 

n = [0:20]; xl = n.*(stepseq(0,0,20)-stepseq(10 , 0 , 20)); 

x2 = 10*exp(-0.3*(n-10)).*(stepseq(10,0,20)-stepseq(20 , 0 , 20)); 

x = xl+x2; 

subplot(2 f 2 f 2) ; stem(n,x); title(' Sequence in Exercise b’) 
xlabel( 'n’ ); ylabel( ’x(n)’ ); 


c. x(n) = cos(0.04rm) + 0.2 w(n), 0 < n < 50. 


n = [0:50]; x = cos(0.04*pi*n)+0.2*randn(size(n)); 
subplot(2 , 2 , 3); stem(n,x); title(’ Sequence in Exercise c ? ) 
xlabel( 'n ! ) ; ylabel( ! x(n) ’ ) ; 


d. x(n) 


5, 4, 3, 2, 1, 5, 4, 3, 2, 1,5, 4, 3, 2, 1, ... 7 , 


-10 < n <9. 


n = [-10:9]; x = [5 A 4 A 3 A 2 A 1]; 

xtilde = x’ * ones(1,4); xtilde = (xtilde(:)) ! ; 

subplot(2,2,4) ; stem (n, xtilde); title( 1 Sequence in Exercise d') 
xlabel(’ n ? ) ; ylabel( 'xtilde(n) ’ ) ; 


Sequence in Exercise a Sequence in Exercise b 




Sequence in Exercise c Sequence in Exercise d 




Figure 1.3 : Sequences in Exercise 2 


Exercise 3: 

Let x(n) = {1,2, 3,4, 5, 6, 7, 6, 5, 4, 3, 2, 1}. 

t 

Determine and plot the following sequences. 

a. xi (n) = 2 x(n - 5) - 3 x(n + 4) 

b. X2(n) = x(3 - n) + x(n) x(n - 2) 

Solution: 

The sequence x(n) is nonzero over -2 < n < 10. Hence: 

»n = -2:10; x = [1:7,6:-1:1 ]; 

Will generate x(n). 

a. xi (n) = 2 x(n - 5) - 3 x(n + 4). 

The first part is obtained by shifting x(n) by 5 and the second part by shifting x(n) by -4. This shifting and the 
addition can be easily done using the sigshift and the sigadd functions. 

n = -2:10; 
x = [1:7,6:—1:1] ; 

[xll,nll] = sigshift(x,n, 5) ; [xl2,nl2] = sigshift(x,n,-4) ; 

[xl,nl] = sigadd(2*xll,nil,-3*xl2,nl2); 

subplot(3,1,1); stem(nll,2*xll);xlabel ('nil') ; ylabel( 1 xl1(nl1) ! ); 
title (' Sequence in Exercise 3. a’) ; axis ([-15 15 -25 20]); 

subplot(3,1,2); stem(nl2,-3*xl2); xlabel( ’nl2 ! ); ylabel( f xl2(nl2)’ );axis 
( [-15 15 -25 2 0]); 

subplot(3,1,3) ; stem(nl,xl) ; xlabel ( 'nl' ) ; ylabel( 'xl(n) ' );axis ([-15 15 
25 20] ) ; 

xlabel( 'n' ) ; ylabel( 'xl(n) ’ ) ; 
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Figure 1.3: Sequences in Exercise 3.a 


b. X2(n) = x(3 - n) + x(n) x(n - 2) 

The first term can be written as x(-(n - 3)). Hence it is obtained by first folding x(n) and then shifting the result by 
3. The second part is a multiplication of x{n) and x{n- 2), both of which have the same length but different support 
(or sample positions). These operations can be easily done using the sigfold and the sigmult functions. 

clear all, clc 
n = -2:10; 
x = [1:7,6:-l:l] ; 

[x21,n21] = sigfold(x,n); [x21,n21] = sigshift(x21,n21,3); 

[x22,n22] = sigshift(x,n,2); [x22,n22] = sigmult(x,n,x22,n22); 

[x2,n2] = sigadd(x21,n21,x22,n22); 

stem(n2,x2); xlabel( 'n2 ' ); ylabel( 'x2(n2) ' ) ;axis ([-15 15 -15 40]); 



















































